Water on Pt(lll): the importance of proton disorder 
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Abstract 

The structure of a water adlayer on Pt(lll) surface is investigated by extensive first principle 
calculations. Only allowing for proton disorder the ground state energy can be found. This results 
from an interplay between water/metal chemical bonding and the hydrogen bonding of the water 
network. The resulting short O-Pt distance accounts for experimental evidences. The novelty of 
these results shed a new light on relevant aspects of water-metal interaction. 

PACS numbers: 68.43.Bc, 73.20.Hb, 68.43.Fg 



The investigation of the interaction of water with solid metal surfaces is of extremely 
high interest in surface chemistry, catalysis and many other technological applications and 
is being carried out extensively by means of experiments and theoretical calculations [l[ Q]. 
In particular, much attention has been devoted to the initial stages of the wetting of tran- 
sition metal surfaces jjj, 4, 5], which occurs through the adsorption of monomers or small 
clusters, and to the structure and stability of very thin water overlayers on Ru(1000) Q], 
Rh(lll) 2|, Pt(lll) B,3 and Cu(llO) 10 ] . In spite of these important efforts, the deter- 
mination of the stable structure of water overlayers is in several cases still controversial, as 
experiments cannot provide conclusive information at the atomistic level. On metals with 
low proton affinity, as in the cases of Pt, Ni and Cu, water arranges in an undissociated 
hexagonal pattern ll|], while on Ru(1000) a half-dissociated monolayer structure has been 
demonstrated to be energetically favored (], 12]. 

The interaction of water monomers and small clusters on Pt(lll) is by far the best 
studied jsl, Q. Water is adsorbed at the top site and the bond with the substrate is 
mainly due to the interaction of the unoccupied 5d states of Pt atoms with the l&i occupied 
molecular state of water. Nevertheless, the periodic adlayer can assume three different 



bilayer arrangements, namely (V39 x V39) J R16.1° (RT39), (V37 x V37).R25.3 (RT37) and 
(y/3xV3)R30° (RT3) , ]l5|, ]l6(, with competing adsorption energies [5j. However ab initio 
calculations showed that the RT39 is the most energetically favored structure 0, [it|] . the 
RT3 has been identified in low energy electron diffraction [18] and STM experiments pjj]. In 
addition, the difference in the adsorption energy between the RT3, RT37 and RT39 structures 
of the wetting layer vanishes to zero as the water coverage increases 5] . Partially dissociated 
OH+H 2 0) structures have only been observed in connection to water production reaction 
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and theoretically studied 
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221 ] . Given this, it is clear that one structure prevails on 



another as a function of the kinetics and the thermodynamic conditions of the experiments. 

While most of the experimental evidences are well described by the remarkable wealth of 
theoretical results in literature, recent XAS and XPS results indicate a water layer distance 
to the substrate (2.4 A) by far shorter than the theoretically predicted one [9|. Indeed, the 
previously reported theoretical studies predict that the most stable water configurations is a 
bilayer, where the shortest Pt-0 distance (2.8 A) is the result of a weak water-platinum bond. 
In this configuration the water network is strongly bound and it is left essentially unchanged 
by the interaction with the metal substrate. This discrepancy between experimental evidence 
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and theoretical results is still an open problem and in this Letter we address this issue via 
a computational approach. We use a large unit cell, since it allows for a less biased search 



of the lowest energy structure as underlined in Ref. 23j, and investigate the possible ro 



of proton disorder which is known to be essential to stabilize not only bulk crystals [33], 
but also crystalline films 24|; moreover it has been already supposed to play an important 
role in determining the wetting structures on transition metals surfaces [3, [^J. We show 
that for this kind of system the proton disorder is essential to understand the ground state 
structure. In addition, for the fully proton disordered system, some of the water molecules 
form a strong bond with the substrate leading to a short bonding distance (2.3 A). The 
formation of these strong bonds leads to a weakening of the hydrogen-bond network and the 
remaining molecules move to a larger distance (3.1-3.5 A). The resulting picture agrees with 
the experimental findings and allows for an interpretation of the adsorption mechanism. 

We use the DFT based finite-electronic temperature method (FEMD) 25] implemented in 
the plane- wave based CPMD code 26] . In this method the electron density and the Hellman- 
Feynman forces are determined via a subspace diagonalization of the high temperature 
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FIG. 1: (Side and top views of the proton disordered H. m ix configuration. The five Pt atoms bonded 
to the oxygens are raised by 0.1 A with respect to the average z-position of the remaining 31 first 
layer Pt atoms. The five bonded oxygen atoms (circled in the top view) stay at 2.33 A from the 
closest Pt atoms, while the rest of the water network lays 3.1-3.4A far from the surface. Hydrogen 
bonds are drown in red (blue) when the 0-0 distance is smaller (greater) than 3.0 A (see text). 
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electron density matrix. The subspace is expanded in a plane-wave basis set that in our 

set-up is cutoff at 60 Ry, a value which, according to our previous work [4|, is sufficient to 

describe satisfactorily the molecule/metal interaction. We use pseudopotentials generated 

I I 

according to Troullier- Mart ins [27| scheme; all the pseudopotentials were accurately tested 
to reproduce bulk and surface properties of the metal, and correct structural properties for 
the molecule. We use the PBE [2^| generalized gradient corrected functional. The prototype 
system consists of a (111) surface of platinum represented by 3 and 4 layers (as explained 
later on regarding the geometry optimization procedure) and a layer of water composed of 
24 molecules. We employ a 6x6 supercell of hexagonal symmetry with the cell dimension 
in the direction perpendicular to the surface equal to 25.5 A so that the thickness of the 
vacuum between the water molecules and the bottom layer of the image slab of Pt(lll) is 
equal to about 15 A; due to the extended size of the box, we performed the calculations 
at the T-point approximation. Geometry optimizations, using the BFGS algorithm, were 
done first with only 3 layers until the maximum component of the ionic forces was below 
2 x 10 -3 atomic units, then we add a fourth layer and further optimize the structure until 
the convergence of the root-mean square force was below 1CT 3 atomic units, the maximum 
component of the ionic force was at least below 3 x 1CT 3 and the energy changes were below 
0.01 eV, i.e. about 10~ 4 eV per molecule into the adsorption energy of water. 

Five different proton arrangements are considered: a fully proton-ordered configuration 
(H^otn)? similar to the ones already studied in literature, two in-plane proton-disordered 
structures (H. up and B.d OW n) and two fully proton disordered structures, one dissociated (H<& ss ) 
and one undissociated (H m j X ). Following the convention of previous works 29| we define 
H up and B. down as the configuration where the hydrogen atoms not participating to hydrogen- 
bonds point toward the vacuum or metal surface, respectively. The initial configuration of 
the water film is the one of an ideal (0001) bilayer of hexagonal ice stretched in the xy plane 
so to match the surface lattice parameter of the metal. In the case of Pt this mismatch is 
8%. The proton disordered configurations have been generated by a Montecarlo method 



301 ]. so to minimize either the global or the in-plane dipole moment. Obviously the z 



component of the dipole moment can be set to zero only if there is no constraint on the 
number of protons pointing toward the metal surface: in this case we have the so called H mix 
film. In the initial configurations the water molecules are on the top sites of the platinum 
surface, the vertical separation between the oxygen atoms in the water bilayer is 0.6 A, and 
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the lower layer is 2.6 A apart from the metal slab. In the partially dissociated configuration 
(R diss ), the dissociated protons sit at the free top sites. In Fig. [I] the final geometry of the 
H-mix system is reported; the geometries of the other systems are those reported in Ref . 3] • 

The total adsorption energies per molecule are computed as E a d s = {E to t — Ep t — 
2AEh 2 o) /24, where E tot is the energy of the optimized Pt+water film system, Ep t the energy 
of the bare Pt surface and Eh 2 o the energy of a water molecule. The values are reported in 
Tab. [J and show very good agreement with those of Ref. as for the H^ n ; moreover the 
H mia , results the most stable. In order to characterize the tendency of the hydrogen bonding 
network to become less tight as a result of the chemical interaction with the surface, we 
monitor the number of "strong" hydrogen bonds. An H-bond is defined as "strong" when 
the 0-0 distance is smaller than 3.0 A. In the initial configurations this distance is 2.98 A. 

In all the undissociated structures, but H mix , the H-bond network stiffens during the 
optimization. All the H-bonds shorten, so to recover the ideal 0-0 distance (~ 2.78 A). 
Since the underlying metal surface act as a topological confinement, this process (i.e. 0-0 
shortening) is achieved via a reduction of the vertical spacing of the bilayer (initially as long 
as 0.9 A). The water bilayers relaxes at a rather large distance from the Pt surface, with 
which it does not interact chemically, as confirmed by the charge density plot in FigJ21 In 
Hdiss both the OH groups and the undissociated H 2 bind chemically to the metal surface 
and form a H-bond pattern where strong (0-0 distance: 2.6 A} and very elongated (3.1 A) 
T-bonds alternate. This picture, confirmed by experiments 3l| and previous calculations 



321 ] . suggest a competition between the chemical O-Pt bonding and the strength of the 



H-bond network. 

This balance determines the adsorption mechanism of the H m j X . adlayer, which is by 
far the most stable of the configurations investigated. Due to the full proton disorder, the 
system is less constrained. The optimal balance is found by sacrificing some of the hydrogen- 
bond energy of the two dimensional water adlayer network (top panel of Fig. [3]) and gaining 
chemical energy with the formation of stronger water-metal bonds. Although the periodicity 
constraints are relaxed, the water molecules maintain the initial ordered arrangement in the 
xy plane, while the distribution of vertical heights is clearly bimodal (bottom panel of Fig. 
[3]). The majority of water molecules (80%) form a flat overlayer where the Pt-0 distance 
is distributed between 3.2 and 3.5 A, while the flat water molecules of the bottom layer 
bind tightly to Pt atoms at a distance of ~2.3 A. The electronic density profiles for the 
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FIG. 2: Electron density profiles for the H m i x case. The four possible configurations of water 
molecules are shown in the left panel (the iso-surface level is 0.270 e/A 3 ). From left to right: a flat 
water molecule in the bottom layer forming a O-Pt bond, a water molecule with a hydrogen atom 
pointing to the vacuum, a water molecule forming a hydrogen bond with the surface and a flat 
water molecule in the top layer displaying no interaction with the surface. Right panel: electron 
density difference maps for one of the flat water molecules in the bottom layer, adsorbed via O-Pt 
bond. The green (yellow) surface indicates depletion (accumulation) of electronic density of 0.036 
e/A 3 . Note that this snapshot compares very well to the monomer adsorption electron density 
difference maps shown in Ref. 

latter molecules (see the right panel in Fig. |2j) shows the formation of chemical bonds with 
the surface which are much stronger than in the other structures. Moreover the electron 
density difference of the Pt-Me bonding for the H m j X . case (left panel in Fig. [2]) shows the 
strong chemical bond to be almost similar to that of an isolated monomer |3j, i.e. the ideal 
conformation for the molecular adsorption. As a chemical bond is established between metal 
and oxygen atoms the hydrogen bonding of the water adlayer is weakened, as proven by the 
elongation of a number of hydrogen bonds occurring during the optimization of the H m j X 
system (Fig. [3]). 

In order to provide further evidence to the interplay between Pt-0 chemical bond and the 
stiffness of the hydrogen bonding network in the B. mix system, we substituted Pt with Ag 
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FIG. 3: Top panels: the number of "strong" hydrogen bonds during the geometry optimization 
process for the H m ix case on Pt (Ag) surface is reported on the left (right) side. Bottom panels: 
Evolution of the distances from the metal surface of the water molecules for the H m j x configuration. 
In the left (right) panel we show the distance from the Pt (Ag) surface of the oxygen of the top 
water layer (dashed line), of the oxygen of the Hdown water molecules in the bottom layer (dotted 
line), and (dashed-dotted line) of the "flat" water molecules, i.e. those circled in Fig. [TJ 



34i | . which has a substantially smaller affinity to the adsorption of water [3j. In this manner 



29] . This means 



the metal-oxygen chemical bond is weakened, and performing the geometry optimization we 
observe the up-shifting of the previously chemisorbed flat H2O molecules and the restoration 
of a tight H-bond network (Fig. [3]). The calculated adsorption energy of —0.46 eV is in 
good agreement with the ones reported for proton ordered systems in Ref. 
that for metal surfaces with low water affinity the issue of proton disorder is not relevant. 

While providing novel insights about the adsorption mechanism, it is also remark- 
able that the H m j X on Pt is the only configuration that accounts for the x-ray absorp- 
tion/photoemission spectroscopy measurements in indicating a relatively strong Pt-0 
bonding, with a layer of oxygen atoms 2.3-2.4 A far from the surface. It must be noticed 
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that, in none of the other undissociated structures considered in this work or in the many 



previously reported calculations |5|, |8|, |9|, 



29| | . such a short Pt-0 bond has been observed. 



In conclusion, we have performed an extensive DFT study of the structure and bonding 
of water adlayers on Pt(lll) surface, where, in particular, we have investigated the role 
of the proton disorder. We have found that on this particular surface the most stable 
adlayer is an undissociated configuration with fully proton disorder arrangement. The higher 
stability of the H m j X configuration is given by the interplay between the Pt-water chemical 
bonding and the capability of the water adlayer to preserve its hydrogen-bond network intact, 
though weakened. We provide the evidence that the main features of the mixed layer/metal 
interaction cannot be extrapolated from the study of homogeneously oriented (all down or 
all up) layers. In this sense, this study sheds new light on the nature of the mechanism 
of interaction of water on metal substrates and thus opens new paths of experimental and 
theoretical investigations into the field. 

L.D.S. acknowledges the RZG of the Max Planck Society for the computational resources, 
L.M.G. the financial support of the von Humboldt Foundation. Useful suggestions from K. 
Kremer are acknowledged. 



E ads O-H-Pt O-Pt H-O-Pt OH-Pt 

R°ltn -0-49 3.1 2.8 

Hdown "0.50 3.2 3.1 

R up -0.48 - 2.9 3.5 

B mix -0.53 3.2 2.3/3.2 3.3/3.5 

ftdiss -0.29 - 2.2 - 2.1 



TABLE I: The adsorption energies (in eV) and the distances between oxygen and platinum (in A) 
are here reported for the 5 systems we studied. The labels of the columns have to be interpreted 
as follows: O-H-Pt the water molecule has a proton that points toward the Pt surface, O-Pt the 
water molecule is flat with respect to the metal surface, H-O-Pt the water molecule has a proton 
that points toward the vacuum, OH-Pt refers to a dissociated OH group. 
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